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Abstract: We discuss several improvements in the detection of at- 
mospheric turbulence using SLOpe Detection And Ranging (SLODAR). 
Frequently, SLODAR observations have shown strong ground-layer 
turbulence, which is beneficial to adaptive optics. We show that current 
methods which neglect atmospheric propagation effects can underestimate 
the strength of high altitude turbulence by up to ~ 30%. We show that 
mirror and dome seeing turbulence can be a significant fraction of measured 
ground-layer turbulence, some cases up to ~ 50%. We also demonstrate 
a novel technique to improve the nominal height resolution, by a factor 
of 3, called Generalized SLODAR. This can be applied when sampling 
high-altitude turbulence, where the nominal height resolution is the poorest, 
or for resolving details in the important ground-layer. 
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1. Introduction 

The success of adaptive optics in astronomy has been demonstrated with 8-10m class tele- 
scopes. Performance has been reported with the 10m Keck II Telescope fl], the 8m Very Large 
Telescope 0, the 8.2m Subaru Telescope O and 8m Gemini North Telescope (4) and oth- 
ers. The performance of adaptive optics systems depends strongly on the characteristics of the 
atmospheric turbulence above the telescope Q. Information about the height distribution of 
the atmospheric turbulence in terms of its strength and velocity can be used to optimize adap- 
tive optic models, and prove the case for future installations. Measurements of the atmospheric 
turbulence can reveal important parameters for adaptive optics, the coherence length, ro, the 
coherence time, To, and the anisoplantanic angle, 9o- Other useful parameters include the outer 
scale, Lq, and the power law, j3, of the power spectrum of spatial phase fluctuations (for Kol- 
mogorov, Lo = °° an d J3 = 1 1 /3). Measurement of these parameters has been emphasized with 
the planned construction of Extremely Large Telescopes (ELT) and new adaptive optic tech- 
nologies, such as Ground Layer Adaptive Optics (GLAO). This has led to numerous campaigns 
to characterize the atmospheric turbulence profile at current or proposed observatory sites, for 
example, the Cerro Tololo campaign [6]. At many sites a significant fraction of the turbulence 
has been found near the ground, which is favorable for GLAO. This is promising as GLAO 
relies on compensating the low altitude turbulence and providing a uniform partial correction 
over relatively wide-field of several arcminutes [7-9]. 

Various techniques are used for turbulence ranging, including direct sensing with microther- 
mal sensors on towers ifTOll or balloons ifTTI . remote-sensing with acoustic scattering (SO- 
DAR) 021, or triangulation of scintillation (SCIDAR) 11131 1 141 or of image motion (SLODAR) 
lfT5l[T6l . These techniques have reached a degree of maturity exhibiting reasonable agreement 
when used together in campaigns ((6), Cerro Tololo campaign). Each technique has its benefits 



and limitations in terms of cost, height resolution, height range, temporal resolution, ease of 
implementation and data reduction complexity. 




Fig. 1. Diagram illustrating the geometry of the SLODAR method for a N=4 system. 6 is 
the double star angular separation. D is the diameter of the telescope pupil and w the width 
of the sub-aperture of the Shack-Hartmann Wavefront Sensor (SHWFS) array. The centers 
of the altitude bins are given by A8h where Sh = w/0. The ground-layer can be analyzed in 
higher-resolution by utilization of double stars having larger 6. The diagram illustrates the 
notation to describe the SLODAR theoretical covariance impulse function, Xi(A, Si, 8 j) as 
well as the notation for Generalized SLODAR (see Sectionf4](. 

The SLODAR technique, Fig.Q] is based on Shack-Hartmann Wavefront Sensor (SHWFS) 
that measures the averaged local wavefront derivative or slope across the telescope pupil using 
an array, (N x AO, of square sub-apertures or lenslets. The 2-D spatial cross-covariance of the 
sub-aperture spot motions (or Z-tilts) of the double star components in each frame are averaged. 
Then a 1-D slice is taken along the double star separation axis and then inverted against the 1- 
D theoretical covariance impulse functions providing an estimate of the C^f(h)dh |fl6l . The 
process of fitting covariance impulse functions allows the estimation of the outer scale, Lq, and 
the power law, /3, of the power spectrum of spatial phase fluctuations. The process of obtaining 
Cj](h)dh information from the observational data is further explained in Section[4] The vertical 
resolution is uniform, given by 8h = w/6 where w is sub-aperture size, or lenslet size mapping 
to the telescope pupil. The highest sampled layer, H max = (N — 1) x 8h ps D / '9 , where N is 
number of sub-apertures across the telescope pupil, with the ground layer denoted as ho = 0. 
The vertical resolution and maximum sample height are scaled by cos(£), where £ is the zenith 
distance. 

The exposure times are typically 4 ms to 8 ms and directly proportional to sub-aperture size, 
related to the wind speed crossing timescales. The sub-aperture sizes are designed to be approx- 



imately equal or less than ro, or ranging from 5 cm to 15 cm depending on the median seeing. 
Sensitivity to higher altitude turbulence is reduced because there are fewer longer baselines in 
the pupil, inability to freeze turbulence due to high wind speeds, turbulence strength is usually 
weak (compared to ground-layer) and covariance impulse response decreases with altitude due 
to propagations effects. This paper will attempt to lessen or mitigate the problems arising from 
these factors, particularly focussing on the propagation effects, removal of mirror and dome 
seeing turbulence and improving the height resolution. 

With a scientific motivation to determine the statistical properties of the height distribution of 
turbulence at medium quality astronomical site, we have been pursuing an extensive SLODAR 
campaign to characterize the atmospheric turbulence at the Siding Spring Observatory (SSO). 
Data taken consists of 7x7 SLODAR instrument on the ANU (Australian National University) 
24" (l"=2.54cm) telescope and 17x17 SLODAR instrument ANU 40" telescope. A portion of 
data taken with the ANU 40" telescope implemented real-time data processing at 200 fps to 
improve the observational data quality and reduce storage requirements (no need to store raw 
camera frames). The ANU 17x17 SLODAR instrument has obtained 8h = 75 m and H lnax = 
1200 m when observing 5 Apodis having separation, 6 = 102.9" with a zenith distance, £ = 
50°. For high altitude sampling, we have obtained 8h = 2400 m and H max = 40800 m when 
observing a Crucis having separation, 9 = 4.1" with a zenith distance, £ = 35°. 

We have discovered a tendency for the usual implementation of SLODAR to underestimate 
the strength of high turbulent layers. This was later confirmed in SLODAR numerical simu- 
lations involving phase screens with Fresnel and Geometrical propagation techniques. To ex- 
amine in further detail, theoretical calculations were implemented for the covariance impulse 
functions Ifl6l using modified phase power spectrum to model Fresnel propagation. We describe 
the propagation effects on SLODAR in Section|2] 

Also discovered was that the majority of the turbulence profiles were dominated by the 
ground-layer or zero altitude contribution, Iiq = 0, found in part to be caused by strong mirror 
and dome seeing turbulence. By applying a high pass filter with cut-off of 1-2 Hz to the tem- 
poral centroid data streams, it was possible to remove the mirror and dome seeing turbulence 
from the ground-layer measurement. However, at this stage we point out that the ground-layer 
at Siding Spring dominates the seeing, particularly so on nights with poor seeing. We describe 
the process of removing dome and mirror seeing turbulence from SLODAR data in Section[3] 

In order to obtain improved vertical resolution, we have modified the instrument to optically 
move the zero height analysis plane from the telescope pupil upwards to fractional heights of 
the nominal height resolution, dh. We call this technique Generalized SLODAR and we report 
on the methodology, numerical simulations and observational results in Section@] 

A full report on the Siding Spring SLODAR campaign, which now covers one week per 
season for 18 months, will be forthcoming in a later publication. 

2. Propagation Effects of High-Altitude Layers 

The retrieval of the turbulence profile was initially described by the method outlined by Wil- 
son lfT5l . as the spatial cross-correlation of the centroid data from star A and star B de-convolved 
with the spatial auto-correlation of the centroid data from star A, see Fig.Q] The global X and 
Y tilts of the double star components, A and B, are subtracted from the centroid data, to remove 
any telescope tracking errors. The method assumes the auto-covariance ("auto-correlation" by 
Wilson |[T5ll ) is the spatial invariant impulse response of a thin layer for all sampling heights. 
This assumption simplifies the data reduction but neglects the effects of the global tilt subtrac- 
tion. Later, Butterley et al. lfl6l included the effects of global tilt subtraction by calculation of 
the theoretical covariance impulse functions, based on the power spectrum of phase fluctuations 
of a thin turbulent layer located at each sampling height. For Kolmogorov turbulence, Butterley 



et al. |[T6l show that global tilt subtraction adds tilt-anisoplanatism which could under estimate 
the strength of the highest sampling altitude by up to 20%. The tilt-anisoplanatism is caused 
by the separation of the projected telescope pupils of star A and B onto the turbulent layer, see 
Fig.ffl 

Butterley et. al. Ifl6l calculations do not take into account propagation effects of the high 
altitude turbulent layers to the telescope pupil where SLODAR analysis is performed. Propaga- 
tion may be included by using a modified input power spectrum of spatial phase fluctuations, 
derived by Roddier flTTl as 

rf(f) = W) cos 2 (nXhf 2 ) (1) 

in which P$(f) is the phase power spectrum with no propagation, in other words as the 
wavefront leaves the layer at height ft. P2(f) is then the power spectrum at the ground, ft = 0. 
The nulls of the propagated power spectrum, P^(f n ) = occurs at spatial frequencies, /„ = 

[(« + 0.5) 1 / 2 ]/ Vy, for integer n = 0, 1, and where r/ = (Aft) 1 ' 2 denotes the layer's Fresnel 
length. This results in a measurable decrease, in the variance of the image motion across sub- 
aperture with size comparable to the Fresnel length. The effect is increased by the removal of 
global tilt in the reduction process, as this eliminates low spatial frequencies. 

Therefore, propagation effects are most important for small sub-apertures such as those em- 
ployed in SLODAR. For example, a turbulent layer at ft = 15 km has a Fresnel length of 8.6 cm 
at a wavelength of 0.5 microns. In our site testing observations at Siding Spring, we used small 
(5.8 cm on the 40" and 8.5 cm on the 24" telescopes) sub-apertures for SLODAR because 
the seeing is often poor, with ro about 8 cm in median seeing. These sub-apertures sizes are 
similar to the 5 cm sub-apertures used by the European Southern Observatory (ESO) portable 
SLODAR system using a 40 cm telescope fl31fl6l . 

In Fig. |2] we compare the effects of propagation for a turbulent layer at if = 7709 m, or pupil 
separation, A = 6, projected onto H, where A = HG/w. The target double star referenced in 
calculations is a Cen with separation, 9 = 9.44". The plots show the corresponding co variance 
impulse function A = 6 calculated with numerical simulations involving 300 phase screens us- 
ing Fresnel propagation (only for propagation effects) and the calculated theoretical covariance 
impulse response using the methodology outlined by Butterley et al. |[T6l with the modified 
power spectrum of phase fluctuations (Eq.Q} (only for propagation effects). The numerical and 
theoretical models for the results in Fig. |2] reference the pupil geometry of the ANU 17x17 
SLODAR instrument on the ANU 40" telescope. It is evident that propagation decreases the 
covariance response peak by ~ 20% causing a broadening effect with increasing height. The 
results of Fig.|2]show an excellent agreement between numerical and theoretical results, hence 
validating the use of Eq.Q~]in theoretical calculations and analysis. 

In Fig. [3] we show the effects of propagation on the theoretical covariance impulse response 
functions for increasing height, H, and for sub-aperture sizes, w — 5.8 cm and w = 1 1 .6 cm. The 
covariance impulse functions are modelled for the ANU 17x17 SLODAR instrument using the 
methodology outlined by Butterley et al. Ifl6l . but with the modified power spectrum, Eq.Q] The 
w = 11.6 cm size sub-apertures are modelled using a telescope pupil with twice the diameter 
(D=2m) compared to the w = 5.8 cm size sub-apertures (D = 1 m), but impulse functions are 
identical as pupil geometry is unchanged lfl6l . The propagation effects in decreasing the peak 
covariance values become more severe for increasing height, H = A8h, where in Fig. [3] have 
the values 5ft = 1.29 km and H max = 20.6 km (A = 16). The propagation effects are lessened 
for the larger sized sub-aperture, w = 11.6 cm, but are still significant. 

In Fig. H] we compare propagation effects on observational data taken 16:19 12 April 2006 
(UTC) with the ANU 17x17 SLODAR instrument on the ANU 40" telescope. The 17x17 SLO- 
DAR instrument has relatively small sub-aperture size, w = 5.8 cm, and therefore susceptible 
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Fig. 2. Comparison of propagation effects on the covariance response function using numer- 
ical simulations (black dots with error bars) using Fresnel propagation and with theoretical 
values (red line with crosses) using the modified power spectrum (Eq.[T). The covariance 
plots are (a) longitudinal no propagation, (b) longitudinal with propagation, (c) transverse 
no propagation and (d) transverse with propagation. The longitudinal direction refers to di- 
rection parallel to double star separation axis, aligned along the x-direction of the SHWFS. 
The transverse direction refers to direction perpendicular to double star separation axis, 
aligned along the y-direction of the SHWFS. The comparisons are for single turbulent with 
A = 6 or height H = 7709 m and normalized by their respective A = or height H = m 
functions (i.e. no fitting involved). Plots (b) and (d) show that propagation effects decrease 
the peak covariance value (Si = 6) by ~ 20% compared to plots (a) and (c). 



(a) (b) 

Fig. 3. Comparison of propagation effects on the normalized theoretical covariance re- 
sponse functions for increasing height (H) and different sub-aperture sizes, w. The co- 
variance plots are (a) longitudinal and (b) transverse. The sub-aperture size w = 5.8 cm 
with no propagation (black line with crosses); w = 11.6 cm with propagation (blue line) 
and w = 5.8 cm with propagation (red line). The plots show that propagation effects are 
lessened by the larger sized sub-aperture, w = 11.6cm, but still significant ~ 30%. The the- 
oretical covariance functions are plotted for even A with H = ASh, where 8h = 1 .29 km and 
Umax = 20.6 km (A = 16). The theoretical covariance functions for the SLODAR model are 
discrete valued, defined for integer values. Si, but plotted as continuous lines for clarity. 
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Fig. 4. A comparison between using the theoretical covariance impulse functions with (a) 
no propagation effects and (b) propagation effects on actual observational data taken 16:19 
12 April 2006 (UTC) with ANU 17x17 SLODAR instrument on the ANU 40" telescope. 
The observational data was taken under excellent seeing, r$ = 18.2 cm and exhibits sig- 
nificant high-altitude turbulence. The inclusion of propagation effects, (b), increases the 
strength of the highest turbulence, H ~ 16 km by ~ 25% relative to (a), in agreement with 
Fig. [3] The high-altitude turbulence H > 15 km causes a steeper increase in the cumulative 
turbulence of (b) by ~ 30% compared to (a) ~ 25%. The high-altitude of (b) appears more 
concentrated than (a). 



to propagation effects in the underestimation of high altitude turbulence. The observational 
data represents an individual run of the double star a Cen, with angular separation, = 9.44". 
The dataset consists of 2000 frames at 200 fps and exposure of 2 ms under excellent seeing, 
ro = 18.2 cm, for Siding Spring. The dataset is selected as an example as it exhibits significant 
high-altitude turbulence, subsequently verified in the temporal-spatial cross-covariance data, 
moving at speeds ~ 20 ms -1 . The centroid data was filtered with a 1 Hz high-pass FIR (Fi- 
nite Impulse Response) filter to remove static mirror and dome seeing contributions from the 
ground-layer A = measurement (see Section Ol.The estimation of the C^{h)dh profile was 
implemented by fitting of the transverse (T) theoretical covariance functions for a Kolmogorov 
turbulence power spectrum, j3 = 1 1/3. From Fig.|4]the inclusion of propagation effects increase 
the strength of the high-altitude turbulence H > 15 km by ~ 25%. 

3. Removal of Mirror and Dome Turbulence 

Mirror and dome turbulence manifests as a distinct separate component in the power spectrum 
of angular tilt. Mirror and dome turbulence appears at low temporal frequencies. The mirror 
and dome turbulence contaminates the ground-layer measurement of the turbulence profile. In 
some cases the mirror and dome turbulence can be falsely over-estimate the contribution of the 
ground-layer relative to the free-atmosphere layers. The effects of mirror and dome seeing on 
observational data are shown in Fig. [5] The observational data is taken with the ANU 17x17 
SLODAR instrument on the ANU 40" telescope. 

The mirror and dome turbulence shows strong correlation in the zero spatial offset of the 
temporal spatial cross-correlation frame data for time lags over 100 ms. Following a suggestion 
by R. W. Wilson (private communication with C. Jenkins |[T8l ) we have found that the mirror 
and dome turbulence can be removed from observational data by applying a high pass filter to 
centroid data with cut-off approximately in the range of 1-2 Hz. 

4. Improving the Vertical Resolution 

The nominal height resolution of SLODAR is 8h = w/ 9, where w is the sub-aperture width 
and is the angular separation of the observed double star. To improve the height resolution, 
8h, assuming fixed 6 and fixed exposure time, T, is to reduce the size of the sub-apertures, 
w. The number of signal photons per T is directly proportional to w 2 so reducing w will cause 
photon starvation and hence restrict observations to the brighter double stars (few in number). 
The minimum useable sizes of w range from w = 5 cm for the portable ESO 8x8 SLODAR 
system 11511 and w = 5.8 cm for our ANU 17x17 SLODAR system. Reducing w will also present 
a number of second order effects, including less sensitivity to high altitude layers due to typical 
high wind speeds and propagation effects. The high wind speeds, v will reduce the variance of 
the observed tilts contributed by increasing the effective sampling distance of wavefront tilts 
from w to Tv when TV > w |fl9l . Propagation effects (see Section |2]i will reduce or null power 
in the phase fluctuation power spectrum at spatial frequencies near the Fresnel length, or near 
w, resulting in reduction in the variance of the observed tilts contributed by the layer. Hence a 
physical limit exists for minimum size of w and therefore a minimum height resolution, 8h, to 
ensure a satisfactory performance of SLODAR. 

Therefore, we propose the concept of Generalized SLODAR that improves the height reso- 
lution of SLODAR by combining measurements taken at regularly spaced SHWFS conjugation 
heights at the nominal resolution, 8h. The SHWFS conjugation heights are a fractional amount 
of the nominal height resolution, 8h. By combining Nq datasets at regularly spaced SHWFS 
conjugation heights, new information is provided about the atmospheric turbulence, and is pos- 
sible to achieve a Generalized SLODAR height resolution, 8h* = 8h/Nc- 
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Fig. 5. Plots (a) and (c) are observational data containing significant amounts of mirror and 
dome seeing that cause over-estimating the contribution of the atmospheric ground-layer. 
Plots (b) and (d) are observational data with mirror and dome seeing removed by appli- 
cation of a high-pass filter with cut-off of 2 Hz to the centroid data streams. Plots (a) and 
(b) are the Spectral Energy Density of the Y-centroid slope data (top) and Y-centroid data 
stream (bottom) for star A sub-aperture index [i = 2, j = 5]. Plots (c) and (d) are AVI ani- 
mations (size: (c) 579 KB and (d) 581 KB) of the spatial-temporal 2-D cross-covariance se- 
quences for offset lags of to 20 frames, each offset lag is 5 ms apart. The cross-covariance 
sequences are for X-centroid and Y-centroid data, and normalized to the zero spatial offset 
peak, [Si = 0, Sj = 0], for zero offset lag, % = ms. The contribution of mirror and dome 
seeing to the ground-layer measurement ([Si = 0, Sj = 0], T = ms) of plot (c) is about 
48%, observed as excess residual for [Si = 0, S j = 0] for T = 50 ms. Plots (a)-(d) reference 
the observational dataset of a Cen consisting of 20 s of data at 200 fps (4000 frames), 
taken 12:43 20 June 2006 UTC with the ANU 17x17 SLODAR instrument on the ANU 
40" telescope. 



We begin this section by describing the process of retrieving an estimated profile of the atmo- 
spheric turbulence from observational data. We then introduce the notation and summarize the 
results of theoretical covariance impulse functions derived in the paper by Butterley et al. lfl6l . 
We then extend their results to the case of Generalized SLODAR by defining a new set of coor- 
dinates and the methodology for combining measurements and retrieval of the super-resolution 
turbulence profile. We validate the methodology by showing the results of a numerical simula- 
tion of resolving two phase screens closely separated in height. We then apply the methodology 
to observational data and demonstrate an improvement by a factor of two in height resolution. 

The SLODAR technique does not measure the atmospheric turbulence profile directly but 
needs to be recovered from the wavefront slope cross-covariance data from stars A and B. 
The atmospheric turbulence profile is an internal property of the system that can be estimated 
by fitting modelled theoretical covariance impulse response functions of thin turbulent layers 
spaced by 8h. This model assigns equal a priori probability to all heights (unbiased model). 
The fitting procedure can be modelled as system of linear equations in matrix form, Ax = b. 
A is the kernel matrix with column vectors corresponding the theoretical covariance impulse 
response functions, b is an ensemble average of the observed atmospheric turbulence covariance 
profile (that also includes systematic and statistical noise), represented as a column vector, x is 
the quantity that we seek, an estimate of the atmospheric turbulence profile, represented as a 
column vector of strengths of each thin layer. 

We must note that the estimate of the atmospheric turbulence, x, is based on the input data, 
b, assumptions made by the model, A, and the process to recover x, (inversion models). The 
system is over-determined as there are more equations than variables so matrix A cannot be 
directly inverted and a least squares solution is sought. The system solution, x, can be found by 
least squares inversions = A + b, where A + is the pseudo inverse of A. However, we note as b 
contains unwanted noise so the system may be un-stable and hence the solution, x, invalid. 

We can improve the inversion model by using the prior information that the layer strengths 
are a positive quantity, x > and that b is possibly corrupted with Gaussian noise. Such an 
inversion model is the Non-Negative Least Squares (NNLS) algorithm. We have found through 
simulation that the NNLS algorithm recovers the input atmosphere model more accurately than 
other regularization algorithms, such as MAXENT and Tikhonov regularization. The NNLS 
algorithm is implemented in MATLAB as the routine Isqnonneg and performs well on compact 
sources (minimal smoothing). Hence the NNLS algorithm suitable with the thin-layer model 
assumption of the atmosphere, as verified with typical high resolution measurements of atmo- 
spheric turbulence ifTTI . 

For SLODAR, the 2-D theoretical covariance impulse response function to a turbulent layer 
at altitude H for the wavefront tilts in longitudinal (L) direction, given that the double star 
separation axis is aligned with SHWFS x-axis, has been derived by Butterley et al. |[T6l : 

X L (A,8i,Sj)=-^— £ C x j i (A) (2) 

*V cross ,. . . . .1 .1 ■■■>> > J 
valid i j, i ,7 

The function d x , , (A) describes the theoretical covariance of x-directional slopes for a 

cross-pair of sub-apertures with lateral pupil spatial offset (8i,8 j) and layer height, H = A8h, 
after global tilt subtraction. The number of cross-pair lenslets having the same lateral pupil 
spatial offset (8i,8j) for a given layer height, H,is denoted by N cr0 ss- The notation used to 
describe the theoretical covariance impulse response function, Xl(A, 8i, 8 j), is shown in Fig.Q] 
The indices [i,j] refer to the lenslet index for star A and [/',/] for star B. A cross-pair of 
lenslets has a lateral pupil spatial offset defined by (8i,8 j)=(i ! — i,f — j), specified in units of 
the sub-aperture width, w. The lenslet index, i, takes on integer values ; = {1,2, ...,N}, where 
N is the number of lenslets mapped across the diameter of the telescope pupil. Likewise for 



indices j, i' and /. The lateral pupil spatial offset, 8i, takes on integer values Si = {1 — N,2 — 
N, ...,0, 1, ...N — 2,N— 1}, specified in units of sub-aperture width, w. Likewise for 8 j. 

If the double star separation axis is aligned with the SHWFS x-axis, then processing is sim- 
plified by considering the covariance function of the tilts in longitudinal (L) or transverse (T) 
directions relative to the lateral pupil spatial offset, (8i,8j), (units of w). A turbulent layer 
height at H corresponds to a lateral pupil spatial separation, A = HO/w of telescope pupils, 
specified in units of the sub-aperture width, w. The lateral pupil spatial separation, A, is an off- 
set of the projected telescope pupils along the x-direction at the layer altitude, H, and takes on 
integer values A = {0, 1,2, ...,N — 1}. The physical separation of a pair of sub-apertures with 
a spatial offset (<5/,<5 j) projected on a layer at height, H, is then (u x , u y ) where u x = \8i + A\w 
and u y — \8j\w, is used by d x . ., ., (A) function. Hence the estimated strengths of the layers 

are defined with height bins of widths 8h and centered at A8h. However, the practical height 
resolution, 8h, may be poorer depending on the signal-to-noise ratio of observational data and 
the inversion model implemented to recover the estimated strengths. 

The 2-D theoretical covariance impulse response function, Xi(A,8i,8j), is an accurate 
model, and takes into consideration the pupil geometry (mapping of circular or square sub- 
apertures on the annular telescope), turbulence power spectrum (/3, L () ) and effects of 'global' 
tilt subtraction (tilt anisoplanatism) required to remove telescope tracking errors. The parame- 
ters A and (Si,8 j) are integer valued and hence Xl(A, 8i, 8 j) is a discrete function that models 
the impulse response of equally spaced thin layers with height, A8h. The discrete impulse re- 
sponse function Xi(A,8i,8 j) is in a format that is compatible with the discrete observational 
covariance profile C^ bs (Si, 8 j). Hence the discrete function Xi(A, Si, S j) can be specified in 
matrix form, A, to model the system as a set of linear equations, Ax = b, and then inverted 
to solve for layer strengths, x = A + b. To further explain the process the discrete observational 
covariance profile C['£ bs (Si, 8 j) can be modelled as a linear equation in the form 

C x L f\8i, 8j) = x X L {0, 8i, Sj) +xiX L (l, 8i, Sj) + ■■■ +x N -iX L (N - 1 , 8i, Sj) (3) 

where, as noted, the combination of 5; and 8 j maps the theoretical covariance impulse re- 
sponse, Xl(A, 8i, 8 j), of a particular height, H = Aw/ 9. Expressing as a set of linear equations, 
Ax = b, where x is a column vector of layer strengths: 



col{X L {0,8i,8j)} col{X L {\,8i,8j)} ... col{X L (N - l,8i,8j)} 



Xff-l, 



,l{clf\8i,8j)} 



(4) 



where co/{x} denotes the process that serializes the 2-D data, x, into a 1-D column vector 
by stacking columns of x with increasing 8i. 

The covariance impulse function is further simplified by taking a 1-D cut along y = of the 
2-D theoretical covariance function, Xi(A, 8i, 8 j), or by setting / = j or 8j = 0. 



X L (A,Si) = -l— £ C x j (A) (5) 

IV cross '■' 

The 1-D theoretical covariance function,Xi(A, 5/) , is calculated for integer valued lateral 
pupil spatial separations, A = {0, 1,2, ...,N— 1} and the condition A = corresponds to com- 
pletely overlapped telescope pupils projected on the SHWFS. This configuration is when the 
SHWFS is conjugated to the telescope pupil (ho=0 km), refer Fig.Q] 

Information between the nominal height bins, Adh, can be found at non-integer lateral pupil 
spatial separations, Q.^ = A + rj^ = {0 + ^, 1 +77^,2 + rj^, ...,N— 1 + Tj/.}, where rj^ takes values 
between and 1, where k is the index of the group of Nq Generalized SLODAR datasets, 
k = {0, 1, ...,A^g — 1}- The value Q.^ can be obtained by moving the SHWFS conjugation height, 
ho, upwards by fractional amounts of the height resolution, = rj^Sh, and is illustrated in 
Fig.Q] Moving the conjugation height, h%, results in a lateral pupil spatial offsets, dm^ = 8i + 
rjk and dm^ = 77^ for 8i = 0, corresponding to a lateral pupil spatial separations, £li = 77^. 
Hence telescope pupils are no longer completely overlapped at the SHWFS but separated by 
a fractional amount of a lenslet. The non-integer lateral pupil spatial separations, Q.^, can be 
thought of sampling new and unique spatial offsets, 8i+ in the telescope pupil. 

The aim of Generalized SLODAR is to reconstruct a super-resolution turbulence profile by 
combining several datasets, k, having unique lateral pupil spatial separations, £1%, and with 
equal height resolutions, 8h. The methodology for Generalized SLODAR is shown in Fig. [6] 
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Fig. 6. The data reduction methodology for Generalized SLODAR 



To provide an unbiased super-resolution profile the lateral pupil spatial separations, £1^, must 



be equally spaced and hence require r\ k to also be equally spaced. The fractional spacings, r\ k 
are then given by r] k = k/N G and therefore rj k = (1 /W G ){0, 1,2, ...,Ng — 1}. 

We denote the observed global-tilt removed covariance profile in the longitudinal (L) direc- 
tion of double star having a separation axis aligned along the x-axis of the SHWFS, for an 
individual Generalized SLODAR dataset, k, to be Cf'f' s (Si). Note that the symbol defined for 

the observed covariance profile C?f bs (8i) should be clearly distinguished from the theoreti- 
cal covariance function for a cross-pair of lenslets, C x . ., ,(A). We now need to transform the 
observed covariance profile from a local lenslet-based coordinate system to a global coordi- 
nate system, C^f s (8m k: ), referenced to rj k = 0, or lateral spatial offsets in the telescope pupil 
at Iiq = 0. The global coordinate system of an individual Generalized SLODAR dataset, k, is 
defined as 8m k = 8i + rj k , specified in units of the sub-aperture width, w. To construct the ob- 
served super-resolution covariance profile, C™(<Sm*), requires the C^' k ° s (8m k ) profiles to 
be first scaled to normalize fluctuations in seeing and then interleaved. The scaling parame- 
ter, a k , normalizes C^ bi (8m k ) to have equal seeing and hence remove any bias effects, and 
defined as 

Ay(&=o) 

h x L f\8i = Q) 

where A^ hs (8i = 0) refers to the peak of the centroid-noise removed auto-covariance func- 
tion for dataset k, and proportional to the total atmospheric seeing. For most cases the scal- 
ing parameter, a k , is close to unity, a k » 1 . The observed super-resolution covariance profile, 
C2' r '° 6s (5m*),isthen 

C^" h \8m*) = [}a k & L f{8m k ) (7) 

k 

where 

8m* = [jSm k (8) 

k 

is the combined spatial offsets relative to the Classical SLODAR, expanded form: 

8m* = {8mo U 8m\ U ... U 8m^ G -i } (9) 
The pupil spatial separations are denoted by: 

Cl k = A + 7] k (10) 

where 

k 

is the combined pupil spatial separations relative to the Classical SLODAR, expanded form: 



Q.* = {£2 U£1iU...UQjv g _i} 
we now need to calculate the super-resolution theoretical function: 



(12) 



xm*,8m*) = ^ I <*y (tf) (13) 

Across ... , / ' 5 5 
vana m.l.m 

where to and Z are now indices that reference a higher sampled SHWFS at fractional spacings 
of a sub-aperture, w* = w/Nq with total samples, ofN* ~NqN. Due to the complexity and time 
required to compute X[ (Q.* , 8m* ) it is best to approximate with interpolation methods. Through 
numerical simulations involving phase screens, it found that cubic interpolation method is suit- 
able for the X[ (£2* , 8m* ) function and spline interpolation for Xj (£2* , 8m* ) function. 

The theoretical covariance function, X£(£l*,5m*), can now be constructed in matrix form, 
A, to model the system as a set of linear equations, Ax = b, and then inverted to solve for 
layer strengths, x = A + b. Due to the larger size of the matrix A, it best to use a positively 
constrained, x > inversion method for compact sources (minimal smoothing to x), such as the 
Non-Negative Least Squares (NNLS) algorithm implemented as the MATLAB iterative routine 
Isqnonneg. 

From theoretical and numerical simulations the technique is successful for Nq = 3. For No = 
6 the results are progressively poorer due to a larger matrix being increasingly sensitive to noise. 

We validate the Generalized SLODAR methodology presented in this paper by showing the 
results of a numerical simulation for Nq = 3 resulting in an effective height resolution, 8h* = 
8h/3. We confirm this by clearly separating two phase screens separated in height by h = 28h*. 



Parameter 


Value 


Description 


e 


9.44" 


double angular star separation 


A 


0.5e-6 


mean wavelength 


D 


1.02 m 


telescope diameter 





0.45 


secondary / primary obstruction ratio 


w 


0.06 m 


sub-aperture width (square) 


N 


17 


number of sub-apertures across telescope diameter 


8x 


1 cm/pixel 


waverfront and pupil sampling 


8x h 


200 m/pixel 


vertical resolution of pupil sampling 


N frames 


4000 


number of independent atmospheric realizations 


Propagation 


Geometrical 


phase screens added together for pupil wavefront 


Hi 


6800 m 


height of phase screen for layer 1 


H 2 


7600 m 


height of phase screen for layer 2 


ro 


0.3 m 


total integrated seeing 


JS 


11/3 


power law of phase power spectrum (Kolmogorov) 


N G 


3 


number of generalized datasets 


m 


0, 1/3, 2/3 


fractional pupil offsets for generalized datasets 


8h 


1200 m 


nominal resolution of each generalized dataset 


8h* 


400 m 


super-resolution of combined profiles 



Table 1 . Parameters for the numerical simulation 



We model the simulation after the double star a Cen and the ANU 17x17 SLODAR instru- 
ment on the ANU 40" telescope. The parameters of the simulation are listed in Tab. Q] The 
Generalized SLODAR is simulated by sequentially moving the layers down in vertical height 
by 8h* or 400m for each fractional generalized pupil offsets, A lateral pixel offset of 1cm 
corresponds to vertical height of 200m. Therefore, for each dataset, decreasing the separation 
of telescope pupils of star A and star B as projected onto the phase screens H\ and H2 by two 
pixels (2x200m) achieved Generalized SLODAR. The wavefronts for each star at the SHWFS 



is calculated by extracting the part of the phase screen that the pupils project on and then adding 
together for each layer H\ and H2. 

The results of the simulation clearly separated the phase screens as illustrated in Fig. [7] 
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Fig. 7. Numerical simulation of Generalized SLODAR using parameters of Tab. Q] The 
objective is to fully separate two thin turbulent layers with height separation, AH = 28h* 
where 8h* = SIi/Nq with 8h = 1200 m and Nq = 3 generalized datasets. Plots (a)-(d) 
denote the numerical results as filled circles with error bars (black). Plot (a) shows the 
averaged longitudinal auto-covariance profile using star A; plot (b) shows the combined 
super-resolution longitudinal cross-covariance profile, C L s (8m*), using star A and star 
B; plot (c) shows the super-resolution C^(h*)dh* profile obtained by fitting super resolu- 
tion kernel, X£(£l*, 8m*), to the cross-covariance profile, plot (b), using 4000 atmospheric 
realizations; plot (d) is similar to plot (c) except using 2000 atmospheric realizations. Plot 
(a) denotes the best theoretical fit with parameters j8 = 3.63 and po = 0.31 ±0.01 m as 
continuous line (red); theoretical G-tilt of sub-aperture as circle (black); theoretical Z-tilt of 
sub-aperture as square (black). Plot (b) denotes the best theoretical fit of super-resolution 
C^(h)dh profile, plot (c), as continuous line (red). Plots (c) and (d) denote the modelled 
atmosphere with parameters = 3.67 and po = 0.3 as stem lines with asterisks (red). 
The numerical simulation results shown in plots (a)-(d) confirm the validity of using the 
methodology outlined in Section|4]and illustrated in Fig.[6] 



We now apply the Generalized SLODAR methodology to observational data by combin- 
ing two SLODAR datasets of the double star a Cen, angular separation of 9.44", at Siding 
Spring Observatory. The first dataset was captured at 10:04 21 June 2006 (UTC) with SHWFS 
conjugation height Iiq = 550 m (rji = 0.5) and second dataset was captured at 10:59 21 June 
2006 (UTC) with SHWFS conjugation height = 990 m (tj 2 = 0.9). The third dataset hav- 
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Fig. 8. Observational data of Generalized SLODAR, Nq = 2 datasets, having fractional 
offsets Tji = 0.5 and 772 = 0.9, and nominal height resolution, Sh = 1 102 m. The double star 
is a Cen, with separation 9 = 9.44", observed 10:04 (77O & 10:59 (772) 21 June 2006 (UTC) 
with the ANU 17x17 SLODAR instrument on the 40" telescope at SSO. The plots denote 
the observation results as filled circles with error bars (black). Plot (a) shows the combined 
super-resolution transverse cross-covariance profile, Cj X '° s (8m*), using 4000 frames for 
each dataset. Plot (b) shows the super-resolution C^(h*)dh* profile obtained by fitting super 
resolution kernel, Xj(Q.* ,8m*), to the cross-covariance profile, shown as continuous line 
(red) in plot (a). The observational data gives the seeing conditions as po = 0.095 ± 0.004 m 
and a power law of /3 = 3. 15 ±0.04. Plot (b) shows a 2x improvement in nominal resolution, 
Sh* ~ 8h/2, for the C^(h*)dh* profile, indicating strongest turbulence is near the ground 
< 50 m. 



ing SHWFS conjugation height ho = m (TJ = 0) was taken 12:43 21 June 2006 (UTC) and 
excluded in the analysis as the atmospheric seeing changed significantly (poor seeing) during 
the lhr45mins of observing downtime. Note the fractional spacings of Tji = 0.5 and 771 = 0.9 
are not regularly spaced but the methodology and results remain valid. Each dataset consists 
of 4000 frames captured at 200 fps using a fixed exposure of 2 ms with centroid sequences 
from each lenslet pre-processed by 1 Hz high pass FIR filter to remove mirror and dome see- 
ing contributions from the ground-layer turbulence measurement bin. The results are shown in 
Fig. [8] and clearly demonstrate an improvement in height resolution by a factor two over the 
nominal resolution of 1 100 m, providing an 'effective' resolution of 550 m. The error bars are 
one standard deviation calculated by dividing the dataset into 10 segments of 400 frames. As 
the SHWFS conjugation height /io = Om (rjo = 0) was excluded from the analysis we added 
a single impulse response function for Q.* = to model the ground-layer. We note that the 
turbulence bin for height 550 m does not register any strength. The error bar for this bin con- 
strains the lowest turbulence to be below ~ 50 m, as otherwise the finite width of the covariance 
impulse response for layers at ~ 50 m would cause spill-over exceeding the error bar. 

5. Conclusions 

SLODAR is a simple and valuable technique, particularly for investigation of the ground layer 
with inexpensive and simple equipment. We have shown that some care needs to be taken in the 
analysis of SLODAR data when small (few cm) sub-apertures are used in the Shack-Hartmann 
wavefront sensor, as Fresnel propagation effects can lead to underestimation of high layers, 
thereby overestimating the importance of the common strong ground layers. We have shown 
that pre-filtering the centroid data stream with a high-pass filter with a cut-off around 1-2 Hz can 



remove mirror and dome seeing providing an accurate atmospheric ground-layer measurement. 
We have also shown that a simple optical technique called Generalized SLODAR can yield 
improved vertical resolution in the ground layer at the same time as measuring the high-altitude 
turbulence. 
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